Direct versus indirect sexual selection: genetic basis of colour, size and recruitment in a wild bird

Indirect and direct models of sexual selection make different predictions regarding the quantitative genetic relationships between sexual ornaments and fitness. Indirect models predict that ornaments should have a high heritability and that strong positive genetic covariance should exist between fitness and the ornament. Direct models, on the other hand, make no such assumptions about the level of genetic variance in fitness and the ornament, and are therefore likely to be more important when environmental sources of variation are large. Here we test these predictions in a wild population of the blue tit (Parus caeruleus), a species in which plumage coloration has been shown to be under sexual selection. Using 3 years of cross-fostering data from over 250 breeding attempts, we partition the covariance between parental coloration and aspects of nestling fitness into a genetic and environmental component. Contrary to indirect models of sexual selection, but in agreement with direct models, we show that variation in coloration is only weakly heritable (h(2) < 0.11), and that two components of offspring fitness-nestling size and fledgling recruitment-are strongly dependent on parental effects, rather than genetic effects. Furthermore, there was no evidence of significant positive genetic covariation between parental colour and offspring traits. Contrary to direct benefit models, however, we find little evidence that variation in colour reliably indicates the level of parental care provided by either males or females. Taken together, these results indicate that the assumptions of indirect models of sexual selection are not supported by the genetic basis of the traits reported on here.

Variation of fibre diameter coefficient of variation and fibre curvature across mohair fleeces: Implications for animal selection, genetic selection and fleece evaluation

The presentwork aimed to determine howthe average fibre diameter coefficient of variation (CVD) and fibre curvature (FC) differences between nine sampling sites vary between sex and flock, to identify differences in variability between sampling sites as a result of between animal and between sire variability and to determine correlations between sampling sites in between animal and between sire variability. Australian Angoras (n = 313) from two farms in southern Australia were sampled at 12 and 18 months of age at nine sites (mid side, belly, brisket, hind flank, hip, hock, mid back, neck, shoulder). Staples were taken prior to shearing at skin level and CVD and FC determined. For each shearing, differences in CVD and FC between sampling sites, how these differences were affected by farm, sex, and sire, and the covariance between sites for sire and individual animal effects were investigated by restricted maximum likelihood (REML) analyses. The median mid side CVD at 12 and 18 months of age ranged from 23.6 to 25.1% but the actual range was 16.8–34.2%. The median mid side FC at 12 and 18 months of age ranged from 14.4 to 18.6◦/mm but the actual range was 10.5–26.3◦/mm. The general pattern for CVDwas for the mid back, hip and neck sites to have similar CVD, the brisket, hind flank and hock sites to have larger CVD and the belly to have smaller CVD than the mid side site. The between animal variation for CVD was lowest at the mid back site. This implies that the mid back would be the most effective site for between animal selection for CVD. Heritabilities for CVD (range at 18 months 0.18–0.30) were only about half the heritabilities for mean fibre diameter in the same study. There was a marked anterior–posterior increase in FC at both farms and with both ages. The results give no clear indication of the best site for between animal selection for FC, other than that the hock should be avoided. Heritabilities for FC are moderate to high (range at 18 months 0.44–0.77) and the genetic correlations are high except for the hock. Thus genetic selection for FC at any site, other than the hock, should be effective for changing FC over the entire fleece. There was more variability between animals than between sites and sires. These results are put into context with associated research on variation in mean fibre diameter and staple length.

Evaluation of mature cow weight: Genetic correlations with traits used in selection indices, correlated responses, and genetic trends in Nelore cattle

Genetic correlations of selection indices and the traits considered in these indices with mature weight (MW) of Nelore females and correlated responses were estimated to determine whether current selection practices will result in an undesired correlated response in MW. Genetic trends for weaning and yearling indices and MW were also estimated. Data from 612,244 Nelore animals born between 1984 and 2010, belonging to different beef cattle evaluation programs from Brazil and Paraguay, were used. The following traits were studied: weaning conformation (WC), weaning precocity (WP), weaning muscling (WM), yearling conformation (YC), yearling precocity (YP), yearling muscling (YM), weaning and yearling indices, BW gain from birth to weaning (BWG), postweaning BW gain (PWG), scrotal circumference (SC), and MW. The variance and covariance components were estimated by Bayesian inference in a multitrait analysis, including all traits in the same analysis, using a nonlinear (threshold) animal model for visual scores and a linear animal model for the other traits. The mean direct heritabilities were 0.21 ± 0.007 (WC), 0.22 ± 0.007 (WP), 0.20 ± 0.007 (WM), 0.43 ± 0.005 (YC), 0.40 ± 0.005 (YP), 0.40 ± 0.005 (YM), 0.17 ± 0.003 (BWG), 0.21 ± 0.004 (PWG), 0.32 ± 0.001 (SC), and 0.44 ± 0.018 (MW). The genetic correlations between MW and weaning and yearling indices were positive and of medium magnitude (0.30 ± 0.01 and 0.31 ± 0.01, respectively). The genetic changes in weaning index, yearling index, and MW, expressed as units of genetic SD per year, were 0.26, 0.27, and 0.01, respectively. The genetic trend for MW was nonsignificant, suggesting no negative correlated response. The selection practice based on the use of sires with high final index giving preference for those better ranked for yearling precocity and muscling than for conformation generates only a minimal correlated response in MW. © 2013 American Society of Animal Science. All rights reserved.

Yellow perch strain evaluation I: Genetic variance of six broodstock populations

As a prelude to strain selection for domestication and future marker assisted selection, genetic variation revealed by microsatellite DNA was evaluated in yellow perch, Perca flavescens, from four wild North American populations collected in 2003-2004 (Maine, New York, North Carolina, and Pennsylvania,), and two captive populations (Michigan and Ohio). For the loci examined, levels of heterozygosity ranged from H-e=0.04 to 0.88, genetic differentiation was highly significant among all population pairs, and effective migration ranged from low (N(e)m=0.3) to high (N(e)m=4.5). Deviation from Hardy-Weinberg equilibrium was regularly observed indicating significant departures from random mating. Instantaneous measures of inbreeding within these populations ranged from near zero to moderate (median F=0.16) and overall inbreeding levels averaged F-IS=0.18. Estimates of genetic diversity, Phi(ST), and genetic distance were highest between Michigan and all other broodstock groups and lowest between New York and Ohio. Genetic differentiation among groups did not correlate with geographic distance. Overall, the patterns of variation exhibited by the captive (Michigan and Ohio) populations were similar to patterns exhibited by the other wild populations, indicating that spawning and management practices to date have not significantly reduced levels of genetic variation. (C) 2007 Elsevier B.V. All rights reserved.

An unbiased genetic screen reveals the polygenic nature of the influenza virus anti-interferon response

UNLABELLED: Influenza A viruses counteract the cellular innate immune response at several steps, including blocking RIG I-dependent activation of interferon (IFN) transcription, interferon (IFN)-dependent upregulation of IFN-stimulated genes (ISGs), and the activity of various ISG products; the multifunctional NS1 protein is responsible for most of these activities. To determine the importance of other viral genes in the interplay between the virus and the host IFN response, we characterized populations and selected mutants of wild-type viruses selected by passage through non-IFN-responsive cells. We reasoned that, by allowing replication to occur in the absence of the selection pressure exerted by IFN, the virus could mutate at positions that would normally be restricted and could thus find new optimal sequence solutions. Deep sequencing of selected virus populations and individual virus mutants indicated that nonsynonymous mutations occurred at many phylogenetically conserved positions in nearly all virus genes. Most individual mutants selected for further characterization induced IFN and ISGs and were unable to counteract the effects of exogenous IFN, yet only one contained a mutation in NS1. The relevance of these mutations for the virus phenotype was verified by reverse genetics. Of note, several virus mutants expressing intact NS1 proteins exhibited alterations in the M1/M2 proteins and accumulated large amounts of deleted genomic RNAs but nonetheless replicated to high titers. This suggests that the overproduction of IFN inducers by these viruses can override NS1-mediated IFN modulation. Altogether, the results suggest that influenza viruses replicating in IFN-competent cells have tuned their complete genomes to evade the cellular innate immune system and that serial replication in non-IFN-responsive cells allows the virus to relax from these constraints and find a new genome consensus within its sequence space.

IMPORTANCE: In natural virus infections, the production of interferons leads to an antiviral state in cells that effectively limits virus replication. The interferon response places considerable selection pressure on viruses, and they have evolved a variety of ways to evade it. Although the influenza virus NS1 protein is a powerful interferon antagonist, the contributions of other viral genes to interferon evasion have not been well characterized. Here, we examined the effects of alleviating the selection pressure exerted by interferon by serially passaging influenza viruses in cells unable to respond to interferon. Viruses that grew to high titers had mutations at many normally conserved positions in nearly all genes and were not restricted to the NS1 gene. Our results demonstrate that influenza viruses have fine-tuned their entire genomes to evade the interferon response, and by removing interferon-mediated constraints, viruses can mutate at genome positions normally restricted by the interferon response.

Environmental selection on transcriptome-derived SNPs in a high gene flow marine fish, the Atlantic herring (Clupea harengus)

High gene flow is considered the norm for most marine organisms and is expected to limit their ability to adapt to local environments. Few studies have directly compared the patterns of differentiation at neutral and selected gene loci in marine organisms. We analysed a transcriptome-derived panel of 281 SNPs in Atlantic herring (Clupea harengus), a highly migratory small pelagic fish, for elucidating neutral and selected genetic variation among populations and to identify candidate genes for environmental adaptation. We analysed 607 individuals from 18 spawning locations in the northeast Atlantic, including two temperature clines (5-12 °C) and two salinity clines (5-35‰). By combining genome scan and landscape genetic analyses, four genetically distinct groups of herring were identified: Baltic Sea, Baltic-North Sea transition area, North Sea/British Isles and North Atlantic; notably, samples exhibited divergent clustering patterns for neutral and selected loci. We found statistically strong evidence for divergent selection at 16 outlier loci on a global scale, and significant correlations with temperature and salinity at nine loci. On regional scales, we identified two outlier loci with parallel patterns across temperature clines and five loci associated with temperature in the North Sea/North Atlantic. Likewise, we found seven replicated outliers, of which five were significantly associated with low salinity across both salinity clines. Our results reveal a complex pattern of varying spatial genetic variation among outlier loci, likely reflecting adaptations to local environments. In addition to disclosing the fine scale of local adaptation in a highly vagile species, our data emphasize the need to preserve functionally important biodiversity.

Evolutionary switch and genetic convergence on rbcL following the evolution of C4 photosynthesis.

Rubisco is responsible for the fixation of CO2 into organic compounds through photosynthesis and thus has a great agronomic importance. It is well established that this enzyme suffers from a slow catalysis, and its low specificity results into photorespiration, which is considered as an energy waste for the plant. However, natural variations exist, and some Rubisco lineages, such as in C4 plants, exhibit higher catalytic efficiencies coupled to lower specificities. These C4 kinetics could have evolved as an adaptation to the higher CO2 concentration present in C4 photosynthetic cells. In this study, using phylogenetic analyses on a large data set of C3 and C4 monocots, we showed that the rbcL gene, which encodes the large subunit of Rubisco, evolved under positive selection in independent C4 lineages. This confirms that selective pressures on Rubisco have been switched in C4 plants by the high CO2 environment prevailing in their photosynthetic cells. Eight rbcL codons evolving under positive selection in C4 clades were involved in parallel changes among the 23 independent monocot C4 lineages included in this study. These amino acids are potentially responsible for the C4 kinetics, and their identification opens new roads for human-directed Rubisco engineering. The introgression of C4-like high-efficiency Rubisco would strongly enhance C3 crop yields in the future CO2-enriched atmosphere.

Pervasive positive selection on duplicated and nonduplicated vertebrate protein coding genes.

A stringent branch-site codon model was used to detect positive selection in vertebrate evolution. We show that the test is robust to the large evolutionary distances involved. Positive selection was detected in 77% of 884 genes studied. Most positive selection concerns a few sites on a single branch of the phylogenetic tree: Between 0.9% and 4.7% of sites are affected by positive selection depending on the branches. No functional category was overrepresented among genes under positive selection. Surprisingly, whole genome duplication had no effect on the prevalence of positive selection, whether the fish-specific genome duplication or the two rounds at the origin of vertebrates. Thus positive selection has not been limited to a few gene classes, or to specific evolutionary events such as duplication, but has been pervasive during vertebrate evolution.

Novel Statistical Methods in Quantitative Genetics : Modeling Genetic Variance for Quantitative Trait Loci Mapping and Genomic Evaluation

This thesis develops and evaluates statistical methods for different types of genetic analyses, including quantitative trait loci (QTL) analysis, genome-wide association study (GWAS), and genomic evaluation. The main contribution of the thesis is to provide novel insights in modeling genetic variance, especially via random effects models. In variance component QTL analysis, a full likelihood model accounting for uncertainty in the identity-by-descent (IBD) matrix was developed. It was found to be able to correctly adjust the bias in genetic variance component estimation and gain power in QTL mapping in terms of precision.  Double hierarchical generalized linear models, and a non-iterative simplified version, were implemented and applied to fit data of an entire genome. These whole genome models were shown to have good performance in both QTL mapping and genomic prediction. A re-analysis of a publicly available GWAS data set identified significant loci in Arabidopsis that control phenotypic variance instead of mean, which validated the idea of variance-controlling genes.  The works in the thesis are accompanied by R packages available online, including a general statistical tool for fitting random effects models (hglm), an efficient generalized ridge regression for high-dimensional data (bigRR), a double-layer mixed model for genomic data analysis (iQTL), a stochastic IBD matrix calculator (MCIBD), a computational interface for QTL mapping (qtl.outbred), and a GWAS analysis tool for mapping variance-controlling loci (vGWAS).

Genetic variation in growth traits and yield of rubber trees (Hevea brasiliensis) growing in the Brazilian state of São Paulo

Analysis of variance and covariance was preformed on growth traits (stem girth, bark thickness, total height gain and rubber yield) of 22 open-pollinated progenies of the rubber tree Hevea brasiliensis from an Asian Hevea collection introduced to Agronomic Institute (Instituto Agronômico, Campinas, São Paulo, Brazil; IAC) in 1952. This progeny trial was replicated at three sites in São Paulo state and it was found that at three years from sowing there was statistically significant variation for girth, bark thickness, height and rubber yield. An individual test sites, values of individual plant heritability for girth ranged from ĥ i 2 = 0.36 to ĥ i 2 = 0.89 whereas values for heritability for progeny means ranged from ĥ i 2 = 0.77 to ĥ i 2 = 0.87. These moderate and high heritabilities suggest that a combination of progeny and within-progeny selection would be effective at increasing girth in this population at individual sites. Across sites, values of individual-plant heritability for girth ranged from ĥ i 2 = 0.36 to ĥ i 2 = 0.47, whereas values for heritability of progeny means girth ranged from ĥ x̄ 2 = 0.77 to ĥ x̄ 2 = 0.87. There were high positive genetic correlations between increased girth and bark thickness suggesting that breeding aimed at increasing girth would also increase bark thickness and possibly height. Copyright by the Brazilian Society of Genetics.

Genome-Wide Analysis Reveals Selection for Important Traits in Domestic Horse Breeds

Intense selective pressures applied over short evolutionary time have resulted in homogeneity within, but substantial variation among, horse breeds. Utilizing this population structure, 744 individuals from 33 breeds, and a 54,000 SNP genotyping array, breed-specific targets of selection were identified using an FST-based statistic calculated in 500-kb windows across the genome. A 5.5-Mb region of ECA18, in which the myostatin (MSTN) gene was centered, contained the highest signature of selection in both the Paint and Quarter Horse. Gene sequencing and histological analysis of gluteal muscle biopsies showed a promoter variant and intronic SNP of MSTN were each significantly associated with higher Type 2B and lower Type 1 muscle fiber proportions in the Quarter Horse, demonstrating a functional consequence of selection at this locus. Signatures of selection on ECA23 in all gaited breeds in the sample led to the identification of a shared, 186-kb haplotype including two doublesex related mab transcription factor genes (DMRT2 and 3). The recent identification of a DMRT3 mutation within this haplotype, which appears necessary for the ability to perform alternative gaits, provides further evidence for selection at this locus. Finally, putative loci for the determination of size were identified in the draft breeds and the Miniature horse on ECA11, as well as when signatures of selection surrounding candidate genes at other loci were examined. This work provides further evidence of the importance of MSTN in racing breeds, provides strong evidence for selection upon gait and size, and illustrates the potential for population-based techniques to find genomic regions driving important phenotypes in the modern horse. © 2013 Petersen et al.

Genetic Diversity in the Modern Horse Illustrated from Genome-Wide SNP Data

Horses were domesticated from the Eurasian steppes 5,000-6,000 years ago. Since then, the use of horses for transportation, warfare, and agriculture, as well as selection for desired traits and fitness, has resulted in diverse populations distributed across the world, many of which have become or are in the process of becoming formally organized into closed, breeding populations (breeds). This report describes the use of a genome-wide set of autosomal SNPs and 814 horses from 36 breeds to provide the first detailed description of equine breed diversity. FST calculations, parsimony, and distance analysis demonstrated relationships among the breeds that largely reflect geographic origins and known breed histories. Low levels of population divergence were observed between breeds that are relatively early on in the process of breed development, and between those with high levels of within-breed diversity, whether due to large population size, ongoing outcrossing, or large within-breed phenotypic diversity. Populations with low within-breed diversity included those which have experienced population bottlenecks, have been under intense selective pressure, or are closed populations with long breed histories. These results provide new insights into the relationships among and the diversity within breeds of horses. In addition these results will facilitate future genome-wide association studies and investigations into genomic targets of selection. © 2013 Petersen et al.

Detecting Loci under Recent Positive Selection in Dairy and Beef Cattle by Combining Different Genome-Wide Scan Methods

As the methodologies available for the detection of positive selection from genomic data vary in terms of assumptions and execution, weak correlations are expected among them. However, if there is any given signal that is consistently supported across different methodologies, it is strong evidence that the locus has been under past selection. In this paper, a straightforward frequentist approach based on the Stouffer Method to combine P-values across different tests for evidence of recent positive selection in common variations, as well as strategies for extracting biological information from the detected signals, were described and applied to high density single nucleotide polymorphism (SNP) data generated from dairy and beef cattle (taurine and indicine). The ancestral Bovinae allele state of over 440,000 SNP is also reported. Using this combination of methods, highly significant (P<3.17×10-7) population-specific sweeps pointing out to candidate genes and pathways that may be involved in beef and dairy production were identified. The most significant signal was found in the Cornichon homolog 3 gene (CNIH3) in Brown Swiss (P = 3.82×10-12), and may be involved in the regulation of pre-ovulatory luteinizing hormone surge. Other putative pathways under selection are the glucolysis/gluconeogenesis, transcription machinery and chemokine/cytokine activity in Angus; calpain-calpastatin system and ribosome biogenesis in Brown Swiss; and gangliosides deposition in milk fat globules in Gyr. The composite method, combined with the strategies applied to retrieve functional information, may be a useful tool for surveying genome-wide selective sweeps and providing insights in to the source of selection.